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Abstract 

A quantum network is an open system consisting of several component 
Markovian input-output subsystems interconnected by boson field chan- 
nels carrying quantum stochastic signals. Generalizing the work of Cheb- 
otarev and Gregoratti, we formulate the model description by prescribing 
a candidate Hamiltonian for the network including details the component 
systems, the field channels, their interconnections, interactions and any 
time delays arising from the geometry of the network. (We show that 
the candidate is a symmetric operator and proceed modulo the proof of 
self-adjointness.) The model is non-Markovian for finite time delays, but 
in the limit where these delays vanish we recover a Markov model and 
thereby deduce the rules for introducing feedback into arbitrary quantum 
networks. The type of feedback considered includes that mediated by the 
use of beam splitters. We are therefore able to give a system-theoretic 
approach to introducing connections between quantum mechanical state- 
based input-output systems, and give a unifying treatment using non- 
commutative fractional linear, or Mobius, transformations. 
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1 Introduction 

The aim of this paper is to contribute to our understanding of quantum feed- 
forward and feedback networks by introducing algebraic rules describing how to 
obtain an effective model for a network starting from the canonical description 
of the component devices as unconnected systems and the prescription of the 
interconnections between these systems. In the open systems approach to quan- 
tum mechanics, a unitary dynamics is given for a quantum mechanical system 
(e.g. atom, optical cavity, quantum dot, etc.) and its environment (e.g. optical 
field). When the auto-correlation time of the environment processes is small 
we can employ quantum stochastic approximations and work with a quantum 
stochastic unitary evolution [18] with associated Heisenberg-Langevin equations 



1 



of motion [T3] . The appropriate way to think of the open system is as an input- 
output system [12] where the input process is a causal field representing the 
environment and the output is the scattered field after interaction with the sys- 
tem. Cascading such systems is a basic example of feedforward [Z][II] More 
generally we can consider feedback connections [H] [3] [5] [13] [E] [IH] [IS] ■ 

The natural generalization of this is to consider a graph with quantum fields 
propagating along the edges and quantum mechanical systems at the vertices. 
The simplest quantum network is a single system with input and output, and for 
Markov models the evolution is described in terms of a Hudson-Parthasarathy 
quantum stochastic unitary adapted process [TB] with the system space as ini- 
tial space and the inputs as noise. There is an alternative description based on 
the Chebotarev-Gregoratti Hamiltonian model [9] [17]. The simplest nontrivial 
quantum network will then be two components cascaded in series, and for finite 
time delays this will be non-Markovian: there are several technical and con- 
ceptual difficulties in modelling this in terms of standard (quantum) stochastic 
calculus. To resolve these issues, we extend the Chebotarev-Gregoratti Hamil- 
tonian model to networks, incorporating the various interconnections and time 
delays into the boundary conditions that define the domain of the Hamiltonian. 
It is the authors' opinion that this offers the only feasible way to address the 
topologically non-trivial situation of quantum feedback induced using beam- 
splitter devices. We will show that the limit of small time delays leads to 
a Chebotarev-Gregoratti Hamiltonian of the type associated with Markovian 
models. From this we deduce the rules for introducing feedback/feedforward 
connections into assemblies of Markovian components starting from the com- 
ponent model description. The natural mathematical language to describe this 
is in terms of non-commutative fractional linear transformations of the type 
introduced by C.L. Siegel. 

The notion of quantum feedback for quantum input-output systems has been 
around in one form or another since the late 1970's and has had a major influence 
on theoretical physics and engineering considerations relating to the rapidly de- 
veloping field of quantum feedback control. Our results give a system-theoretic 
approach to introducing feedback. An important step towards a general theory 
of feedforward and feedback connections was made in the papers of Gardiner 
|11| and Carmichael [7] who considered quantum optical networks consisting 
of cascade-connected components with no gauge couplings, and Yanagisawa 
and Kimura [22], [23j who studied the situation where the plants are multi- 
dimensional oscillator systems and the external inputs are Bose fields coupling 
to the plants via emission/absorption interactions. Yanagisawa and Kimura 
were able to exploit the linearity of the dynamics and apply transfer function 
techniques to the resulting networks. The present paper deals with general 
quantum dynamical systems with gauge couplings and is not restricted to linear 
systems. The field channels are assumed to carry quantum stochastic signals 
that satisfy the canonical Ito table. 

In section 2 we review the quantum stochastic and Hamiltonian models 
for open quantum systems, and provide some structural results concerning the 
parameters used to define the models. Our general Hamiltonian description of 
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quantum networks is given in section 3. In section 4 we show how edges can 
be eliminated to provide simpler Markovian network models. Key examples are 
given in section 5, and concrete topological rules are inferred in section 6. 

2 Quantum Markov Input-Output Components 

The concept of quantum input-output systems originates from two independent 
sources: the quantum theory of filtering where the output field is the object 
of indirect nondemolition measurement [2] [3] , and the theory of quantum am- 
plifiers [8]. The latter theory was generalized to a quantum network by Yurke 
and Denker [25j where the network is second quantized quantum wire model 
with quanta satisfying a Klein-Gordon equation with Kirchhoff boundary con- 
ditions at the vertices (though there is no quantum system associated with 
the vertices!). Similarly, starting from a Lagrangian formulation, Gardiner and 
Collett, cf. [H], developed the theory of quantum electromagnetic input-output 
fields interacting with a quantum mechanical system at the origin. We may con- 
sider this set up as the simplest network consisting of an input channel and an 
output channel meeting at a vertex (the system). In the dipole approximation 
they derive a Langevin equations for the canonical observables of the system. 
Following a rotating wave approximation, and a low frequency limit they obtain 
a quantum white noise theory suitable for quantum optics models and which 
is formally equivalent to a Hudson-Parthasarathy quantum stochastic evolution 
for a quantum diffusion |18| . The formulation of the boundary conditions at the 
vertex is then a crucial aspect of the model prescription and relates the output 
field to the input and the system degrees of freedom. 

2.1 Quantum Stochastic Process Description 

The system is modelled as quantum mechanical with Hilbcrt space f). The inputs 
and outputs are carried along semi-infinite quantum channels and are modelled 
as quantum field signals. Specifically, we model the field quanta as propagating 
along the channels with constant velocity c in the direction specified by the 
arrows. We may then parameterize both the input and output lines by a single 
geometric parameter t measuring the arc-time taken to reach the system. A 
single component is sketched below as a two port device having an input and 
an output port. 
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Figure 1: input-output component 

The inputs correspond to the half line IR+ = (0, oo) as they have yet to reach 
the system {t = 0), while the outputs correspond to Mr = (— oo, 0) as they have 
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already passed through the system. Signals therefore have state space L| (R, dt) 
where A is a fixed Hilbert space called the multiplicity space. We shall generally 
consider ^ = C" which means that we have n distinguishable particles. We shall 
consider an indefinite number of these quanta in the wire so that the Hilbert 
space is the Fock space 

5^r{Ll (M, dt)) , 

where F (•) is the bosonic Fock space functor. Note that Fock spaces have the 
functorial property 

r(i)i©f)2) = r(i)i)§5F(i)2), 

and so 5 = ^in ® S'out where 5^i„ = F (L| (M+, dt)) and ^out - T (R-Jt)). 



2.2 Process Description of a Single Markov Component 

Let us fix the multiplicity number as n and take {cj : j — ■ ■ ■ ,n} as a basis 
for R = C". We denote by A, {t) A A{e^(» l[o,t]) and A] (t) = At (e^ l[o,t]) 
the operators describing the annihilation and creation of a quantum in the ith 
channel over the time interval [0,t], respectively. The operator describing the 
scattering from the jth channel to the jth channel over the time interval [0, t] is 
denoted by Ay (t). In particular, Ni (t) = {t) is the observable corresponding 
to the number of quanta in the zth channel over this time. 

We now consider a quantum stochastic evolution as a unitary adapted process 
{V (t) : t > 0} on t) ^ in the sense of Hudson and Parthasarathy [T8] arising 
as the solution to the quantum stochastic differential equation (QSDE) 

dV = idG)V, 1/(0) = 1, 

where 

dG{t) = {S,j - S,j) dAij {t) + L, dA] (t) 

-L^Sy «) dAj (t) - (ilJU + M) ^ dt (1) 

With 5= Sij (g) \ei){ej\ unitary on B (f) (g) j^), i= U (g) (e^l G S ((] (g) M., ()),and 
H G i? (f)) sclf-adjoint. It is convenient to write these as matrices 

Li 

We introduce the processes 

jt{x) ^ vHt)[X(g>i]v{t), 
A,it) ^ vUt)[i(^A{t)]vit). 
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for X e B (f)). They satisfy the QSDEs 

djtiX) = jt{C{X))dt+jt{S]^[X,L,])dAl+jt{[LlX]S,,)dA, 

+jt(sL(X-l)Sfe,)rfA,,(t), (2) 
dA,{t) = jt{S^j)dA,it)+jtiU)dt. (3) 

where C (X) = [X, U] + i x] U - i [X, H]. 
2.3 System Parameters 

The triple (S, L, H), which determines the model, is referred to as the set of sys- 
tem parameters. The coefficients of the QSDE are assembled into the following 
square matrices of (1 + n) dimensions having entries that are operators on [): 



iLtL-iH -LtS 
" L S-1 



(4) 



V - G+n^(-^;L-^H (5) 

M ^ i + nG = i- n + nv=(^|^ , (6) 

where n = ( ^ ? ) ■ We refer to G as the ltd generator matrix of the unitary 
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evolution, and V as the model matrix. The matrix M is called the Galilean 
transformation associated with G. 

Definition 1 Let t) and A be fixed Hilbert spaces. The classes of ltd gener- 
ator matrices 0(f), .ft) and model matrices 9Jt((),.ft) are collections of opera- 
tors G,V G S (f) (g) (C .ft)) of the form ^ and (0) respectively with respect 
to the decomposition f) (g) (C ® .ft) = f) ® ([} .ft) for S € B ([) (g) .ft) unitary, 
L e B (t), t) (g) R), and H e S ([)) self-adjoint. 

It is convenient to set ^1"° = t, = ^] , A^^ = Ai and A'^ = and write 
dX — XaijdA"-^ for a general stochastic integral. We adopt the convention that 
repeated Greek indices are summed over 0,1, - • • , n. The coefficients Xap can be 
assembled into a matrix X of adapted entries. We may compress the quantum 

Its table down to dA^'^dA^" = S^'^dA"'^ where S"'^ is the Hudson-Evans delta 
which equals unity when a = (5 Cz {1, • • • , n} and vanishes otherwise. Note that 

are just the coefficients of the matrix 11. Given stochastic integrals X, Y 
with matrices X, Y respectively, the coefficients of the product XY then form 
the matrix XY + XY -f XHY. 

We may write the generator of the stochastic evolution as dG (t) — Ga/jdA"'' (t) 
and the isometry and co-isometry conditions are 

G + Gt + GtnG = = G + G^ + GnGt. (7) 
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The expression for G above then gives the general solution to this equation. 
The Heisenberg equation for an initial operator X G i? (f)) is then 

J, (X) = X 1 + / js {L^p (X)) ® dA'^" (s) 

and jt {Cap {^}) are the components of a matrix jt (XG + G'^X + G^IIXIIG). 
The super-operators Caf3 are known as the Evans-Hudson maps. 

The output processes are then (t) :— V (f)^ A"'^ (t) V (t) and we deduce 
that 

dA"^ it) ^ jt {Mi^Mfj,) dA^-^ (t). 

This invariance of time, dt = dt, is the motivation for the term "Galilean trans- 
formation", and we now explore some of its properties. 

Definition 2 The Galilean group 6al([),.^) is the group of operators of the 
form (0) m B (t) (g) (C e M.)) where S G i? (f) (g) J^) unitary, L e S ((],() (g) J^) . 

The group identity is I and we readily observe the group laws 

/ 1 \ / 1 \ 

V L2 S2 J L1+S1L2 S1S2 J ' 

[ I S ) - [ -StL St J ■ 
Proposition 3 //M e ©a[(t), j?) then we have the identity MnM^ = IT. 

Proposition 4 The set © is invariant under the action G 1-^ N^GN for 
allN£ ©al(f),j?). 

Proof. Let G e 25 (f),J?) and N e ©a[(f),J?). Setting G' = N+GN we see 
that 

G' + G't = Nt (G + Gt) N = -N+GHG^N 

= N^GNnN^G+N = G'nG't 

since NHN^ = n. Similarly, G' + G't = -G'tRG'. ■ 

Proposition 5 Let G G ©(t),^) with associated Galilean transformation M 
= 1 +nG. Then 

M^GM = G. 
Proof. This follows from the observation 

(I + G^n) G (1 + nG) = G + (Gt + G + G^nG) nG. 
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2.4 Hamiltonian Description of a Single Markov Compo- 
nent 



Let us consider the strongly continuous one-parameter unitary group Uq (t) per- 
forming the time shifts. For instance, taking e (/) to be the exponential vector 
with test function /, we have the action 

Uo{t)u(E)e{f) ^u(E)e{^-tf) 

where dtf (•) — / (• — t). The family {V (t) : t > 0} the forms a right unitary co- 
cycle with respect to Uq, that is V {t + s) = Uq (-s) V {t) Uq (s) V (s) {t, s > 0), 
and we obtain a strongly continuous unitary group U by setting 



r UQ{t)v{t), t>o, 

" 1 i-t) Uq (t), t< 0. 



The generators of U and Uq will be denoted as H and Hq respectively. The 
problem of characterizing H for the class of Hudson and Parthasarathy quantum 
stochastic evolutions has been carried out only in relatively recent times [9] [17] 
and we now recall its explicit construction. 

A vector ^ £ will be a sequence {^m)^^Q where $,„ = $„i (ii, • • • , tm) is a 
f)(g) ^'"-valued function completely symmetric under interchange of its arguments 
and such that 



y- 



\'^in {tl, - ■■ ,tm)\\l^s^ dti ■ ■ -dtjn < OO. 



Let = M\ {0} = (-00, 0) U (0, oo) and define the domain W (M*, J^, f)) con- 
sisting of vectors such that each is differentiable in each of its arguments 
and that 



i) Er=i € ® ^1 (^'"), for each m, 

iii) limt„+i-o± Em 7^ /r™ II*™ (^ir • ' , ^m, Wi)ll^»^-+i exists 



2 

dti ■ ■ ■ dt,n < oo, 



(Note that the left and right limits above need not coincide!) 

The following operators Hq, a {t) for t 7^ 0, and a (0^), are introduced on the 
domain 



.+ 1 



where j is the trace operation from () (g) .^™+^ down to f) (g) ^™ with respect to 
Ci € see [17] for more details. The operator Hq on the Sobolev-Fock domain 
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W (R*, ^, f)) coincides with the generator Hq of translation by time shift on the 
dense subset for which the right and left hand limits in iii) agree. More generally 
we have the relation, which is a consequence of integration by parts with a jump 
discontinuity at the origin, for $, G (K*, .ft, f)) 

n n 

imo^) = (i?o<fi*) (0^) (0") *) -*E(«j- (0^) (0^) *) 

with the sum over an arbitrary orthonormal basis {ci} for the multiplicity space 
Note that we may formally write [3] 




Next, fix a subset D^ c. (S, L) of W (R, ,.ft, f)) consisting of those vectors $ 
satisfying the boundary conditions 

a,(0-)<&-S,feafc(0+)$ + L,$. (9) 

The boundary condition can be written as a Galilean transformation 

a„ (0-) $ = M„^a^ (0+) $ 

where we include the time case Uq — 1. 

Theorem 6 (Gregoratti) f_Z7| / T/ie Hamiltonian K associated with the quan- 
tum stochastic unitary process having the parameters (S, L, H) has dam (H) n 
W (R,, .ft, f)) — Db.c. (S, L) and here it is given by 

H<i>=(^Ho + H- '-L]L, - iL]Sjkak (0+)^ $. 

The Hamiltonian H is essentially self-adjoint on this domain. 

For our purposes, it is most convenient to write the equation for H and the 
boundary condition in terms of the model matrix V as 

- iH^ - (Voo + Vofeflfc (0+) - iHo) $ = {yopap (0+) - iHo) 
a, (0-)$ = V,o$ + Vj-fefl (0+) $ = Vj/3a/3 (0+) $. (10) 

3 Quantum Networks 

The situation sketched in figure 1 will be our simplest example of a quantum 
network: a single component system with one input channel terminating at an 
input port and one output channel starting at an output. A general quantum 
network will consist of several such components connected together and will 
typically have time delays, feedforward and feedback connections. The descrip- 
tion will exhibit separate algebraic, topological and geometric content. While 
the abstract definition is rather involved, figure 2 below gives an example of the 
class of configurations that we wish to consider 
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figure 2: A quantum network 

We first list the basic features of a network. The network consists of a 
collection of components C that arc interconnected in a manner to be described 
shortly. Every component C will have at least one input port and one output 
port. Let us write Vm (C) and Pout (C) as the set of input and output ports, 
respectively, for component C. For each r G Vm (C) we have an associated 
space which is the multiplicity space of the incoming channel. With a 
similar notation for the output ports, we impose the constraint 



reViniC) 



K 



which means that the total multiplicity space of all inputs into a component 
equals the corresponding output one. The sets of all input and output ports in 
the network are then Vi^ = Uc'Pin (C) and Pout = UcPout {C) and we also have 







(11) 



Let riin (C) and riout (C) be the number of input and output ports respec- 
tively in a component and set riin = X)c'^in(^)' '^out ~ X)c^out(C')- Note 
that riin (C) and riout (C) are both non-zero though they need not coincide - 
though the total input/output multiplicities (dim.^c) must be equal! External 
fields propagate into the network along input channels terminating at some of 
the input ports. Likewise, the output fields propagate along output channels 
starting at some of the output ports. Internally, we also have pairs of input and 
output ports connected by further channels. In this way every port is connected 
to exactly one channel. The set of interconnections is described by fixing subsets 
7?.in C Pin and 7?.out C Pout of equal size and a bijection a : 7?.out P-in- The 
pair (s,(t(s)) then determines an internal channel from output port s to input 
port r = a{s). 

Topologically we think of channels as edges and will frequently refer to them 
as such, and distinguish the input, output and internal edges. Each internal 
edge can be written as a pair e — {s,r) - s is the source and r is the range - 
and we will have a corresponding fixed multiplicity space ^ associated with the 
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channel, and this must agree with both ^^^^ and providing an additional 
constraint. The remaining ports are Qm = Pin\7^in and Qout = ^out\^out- For 
each r G Qinwe have a semi-infinite input edge terminating at that input port, 
and similar for each s G Qin- We note the identity 

re Qin seQout 

which implies that the total multiplicity of all external inputs equals that of all 
the outputs. The complete set of edges, including both internal and external, 
will be denoted as £. 

(In the network sketched in figure 2, we have three components C = {Ci, C2, C3} 

with Pin (Ci) = {ri,r2}, Pout (Ci) = {S1,S2}, Pin (^2) = {rg}. Pout (^l) = {^3} 

and Pin(C3) = {rijrs}, Pout ((^3 ) = {s4,S5}- The total network has Pin = 

{ri, r2, ra, r4, rs}, Pout = {51.52,83,54}. The interconnections are via internal 
edges ei = (si,ri) and 62 = (s4,?'2) and so 7?.in = {ri,r2}, 7?.out = {si,S4} with 
o-(si) = ri and cr(s4) = r-2, and Qin = {r3,r4,r5}, Qout = {s2,S3}-) 

In addition to this topological description, we also provide the arc-time taken 
to travel along each internal channel: this determines the various time delays in 
the network. We geometrize the edges by applying a local arctime coordinate 
to each one. Let e = (s,r) be an internal edge then we may assign an arctime 
parameter tg with range {Tj.,Ts) with Ts — Tr being the time taken to travel 
from output port s to input port r. For input channels terminating at r e Qin 
we have an arctime parameter over the semi- infinite range (T,., oc) and likewise 
for output channels leaving s G Qout we have an arctime parameter over the 
semi-infinite range (— 00, Tg). The one-particle Hilbert space for the field quanta 
in an edge e e f of the network is then 

L\r {Tr, 00) , e an incoming edge terminating at r G Qin, 
L|s (— oo,Ts) , e an outgoing edge starting at s € Qout) 
-^^e (^'•(e)) Ts(e)) ! e = (s (e) , r (e)) an internal edge. 

We set 

L'iS)^^L%^ (e) 

The Fock space over all these spaces will be denoted as 5f and by the functorial 
property factors as 

The Hilbert space for the entire network M will then take the form 

^N = ^® d£, 

where i) is the Hilbert space for all the quantum mechanical degrees of freedom 
of the network components. 

In figure 2, we have only sketched the interactions between components that 
are mediated by the channels. Even though it may appear that the network is 
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disconnected, the components may still be coupled, say by a Hamiltonian inter- 
action. We have also stopped short of requiring that the space of components 
factors as, say, f) = (8)cf)C and it convenient not to impose this at this stage. 

Up to this point, we have described the flow through the channels from 
output to input ports. It still remains to describe the trans-component flow. 
This involves the boundary conditions relating the inputs to the outputs at each 
component. The most convenient way to describe this is through the notion 
of the model matrix. With each component C we associate a model matrix 

Vce9Jl(f],ilc). 

Definition 7 Let V,; G *H(f),.^i) for i — 1,2, then the concatenation of model 
matrices is VifflV2 G 2Jt(f),.fti © ^2) defined by 

V Ll Si y'™^, L2 S2 

-iLlLi-iL|L2-z(Hi + H: 
Ll 
L2 

Given the set of component model matrices we may then define the network 
model matrix to be 

net } • 

The network model matrix V relates input ports to output ports and is in- 
troduced independently of the interconnections which go from output ports to 
input ports. It takes the standard form 

/ _iLtL-^H -LtS 
L S 

and, explicitly, for r e Vm and s e Vout, we have components 



ll tl 

2'-2'-2 - 


lH2 - 


L2 




-LlSi 


-L^S2 ' 


Si 








S2 



s 6 Pout 

1 

2 



Voo = -J E 



seVa 



It is convenient to adopt the following block matrix representation with respect 
to these decompositions pT|) of .^nct 










^ L„„„, ) 
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We require that the operators Ssr are contractions with S unitary on f) (Ei^totai, 
and so 

It is possible to use the definition of concatenation in reverse in order to 
analyze a network model matrix into irreducible components. This process is 
substantially more complicated as we have to consider all decompositions (jlip . 

A vector ^ E S^^f can be represented as follows: for each integer to > we 
choose TO locations on the edges £ with local arctime coordinates {ti, • • • , tm}, 
say with tj on edge e (j), to get a vector $,„ {ti, - • ■ , tm) G {®™=i^e{j)) which 
we view as a subset of f) (g) (^'"), where £™ denotes the TO-fold Cartesian 
product of £. We now generalize the class of Sobolev-Fock vectors to networks. 

Definition 8 The class W (M) of Sobolev-Fock functions over a network M is 
the set of vectors ($m)^=o such that each is differentiahle it each of 

its arguments and that 

i) Er=l € t) ® L2 (^^)^ j^^ ^^^^ 



< oo, 



For each internal edge e = (r, s) we introduce the local annihilator density 
Oej (i) on the domain of Sobolev-Fock associated with the edge and the jth 
vector of an orthonormal basis {e^} for the multiplicity space .Se-and < 
t < Ts- We also consider the one-sided limits r+ and and collect into the 
column-block vector of annihilators 



ae{t) = 



ar{T+) = ae(T+), as{T^ ) = a^iT^ ), 



where n — dim.fte. We then define the operator for each internal edge e — (r, s) 
Hoie)=J2j ' {t) i^^a^,, {t) dt = J \l {t) i^^a, {t) dt. 

For the semi-finite external edges we have the corresponding expression with 
either upper or lower limit extended to infinity as appropriate. The total con- 
tribution is then Hq — Eesf (e) and this generalizes the operator ([8]). We 
will denote the ampliation of these operators up to f) (8) 5 by the same symbol, 
and recall our earlier convention that oq = 1. 
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Definition 9 The Hamiltonian for a quantum network M with model ma- 
trix V is the operator on Sjj\j- given by the specification that all Sobolev-Fock 
vectors $ in its domain satisfy the system of boundary conditions 

a, (T-) $ = V..a. (T+) 

re{0}UVi„ 

for each component s € Vout, and on such vectors we have the action 

re{0}U-Pi„ 

Lemma 10 The operator Hj^ in the definition of a quantum network Hamilto- 
nian is symmetric on the restricted domain dam (Hj^) n W (TV) . 

Proof. We now have the network integration-by-parts formula 
(*|iro*)-(iro$|*) = ^ Y {as{T-)^\as{T,-)^)-tY{ar{T+)<^\ar{T+)^) 

for Sobolev-Fock vectors Hj^ is symmetric on the set of Sobolev-Fock 

vectors satisfying the boundary conditions, indeed, 

($|i?Ar*)-(i?A^$|*)=i($|Voo*)-i(Voo$|*)+z J2 (as {T-) <i>\as {T,-) ^b) 
+i Y {(*|Vora. (T+) *) - (Vo,a, (T+) $|*) - (a, (T+) $|a, (T+) *)} 

and, substituting in for Og (T's ) ^ and (Tg^) $ and using the identities Voo + 

Vjo + Y^seVo^i. ^so^so = 0, Vor + Y^seVo^t ^lo^sr, we find the right hand side 
vanishes. ■ 

We may rewrite these equations alternatively as 

a, {T-) $ = ^ S,rar (r+) $ + L,$. 

We wish to identify H_f^ as the physical Hamiltonian for the network, that is, 
show that it is essentially sclf-adjoint. This is a reasonable assumption given 
that it is a direct second quantization of quantum graph models encountered in 
the literature. At the moment, we do not have a proof of this fact but formulate 
it as a conjecture. 

Conjecture 11 We now assume that defines an essentially self-adjoint 
operator. 
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The Hamiltonian Hf^ contains all the physical information about the net- 
work, including the interactions and interconnections. Generally speaking, the 
wave operator V = UqU will not determine a quantum stochastic evolution 
of Hudson Parthasarathy type as the model is typically no longer Markovian. 
We shall show that, in the limit in which time delays along the internal paths 
vanish, we recover a Markovian model which is easily identifiable from H and 
which will have the Chebotarev-Gregoratti form. 

4 Eliminating Internal Edges in the Zero Time 
Delay Limit 

In this section we show how simpler Markovian models can be obtained by 
eliminating edges in a zero time delay limit. We achieve this by showing first how 
to eliminate one edge, and then showing that all edges may be eliminated, with 
the final Markovian model independent of the order in which the eliminations 
were performed. 

Theorem 12 Let eo = (?'OjSo) he an internal channel with time delay To = 
— Tr„ > i] in a quantum network M for which 1 — \/soro *s invertible. In the 
limit To ^ 0^, the network reduces to Mred in which the input and output ports 
are VinX {I'o} and Vout\ {^o} ond the edge eo eliminated. (In the case where tq 
and So are initially in different components, then the components merge.) The 
reduced model matrix V'"^'' then has the components 

K}3 = + y^ro (1 - V.o.o)-' V,„0, (12) 

for /? e {0} U r„,\ {ro} and a G {0} U Vout\ {so}- 

Proof. For $ in the class of Sobolev-Fock vectors we have that a^g $ 
will be norm convergent to a^o (T'^j) ^ as tq ^ O"*". In this limit we therefore 
identify the values a^o (^^t^ ) ^ and (T~^) $ at the start and end of the edge be- 
ing eliminated. The boundary condition Ago (-^sq) ^ ~ ^--'0&{o}yj'Pi 
may then be rewritten as 

ar, (r+) i. = (1 - v,„,j-i Yl v,„.a, (r+) $ 

r£{0}UVin\{ro} 

and substituting into the boundary condition for s € Vout \ {so} yields 

re-Pi„\{ro} 

E v-^, (r+) $ 

re{0}U-Pi„\{ro} 
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The reduced Hamiltonian is then defined by 



Vo,.a. (r+) $ - zi/o-'i* 

re{0}UPi„ 

J2 Vo.a. {T+) $ + Vo.„a.„ (7:+) $ - 

re{0}U-Pi„\{ro} 

J2 v-V(r+)$-.i/-^$, 

re{0}U-Pi„\{rn} 



where TJ^^^ 



Eee£\{eo}-^0(s)- 

The above argument rests on the fact that \\arg (T^) $ — a^^ (r^^) $|| ^ 
as To for Sobolev-Fock vectors $. More expUcitly, by rescahng the 
local time-coordinate on the edge to have fixed range [0,1], we encounter a 
sequence of Hamiltonians H = H'^'^°'> with common domain in Let J7^'^°^ 
and U'^'^'^ be the one-parameter unitary groups generated by the i?'^'^") and 
^red respectively. The Hamiltonian H'^'^'^ has domain in S^j^t^a C with 
the embedding implemented in the obvious way by taking zero quanta in the 
eliminated edge. For each $ £ dom (i7''°'^) we may construct a sequence $("^0) g 
dom (/ff"^")) strongly convergent to $ such that 7j('^o)$(-ro) converges strongly 
to iJ''®'^^. By the Trotter-Kato theorem, e.g. |TD] theorem 3.17, we then have 
that 



lim sup 

To^O o<t<T 



(ro) 



UT 



= 0, 



for all T > and $ e jo^,cd . 



The reduced model matrix is closely related to fractional linear trans- 
formations, see for instance [M]) [IS]) [23 Chapter 10], and we introduce appro- 
priate definitions in the present context. 

Definition 13 Let be the multiplicity space of an edge e = (s, r) so that 
is a subspace of the network multiplicity ^net md let X € B {^e)- The feedback 
reduction o/V e VJl{t),A) through the edge e, with gain X, is the map 

: {Y,X)^T,{V,X) 
where in terms of block decomposition we have 

(V, ^ y^fi + V„,X (1 - \/,rX)-' y,f3, (13) 

where the indices are /3 G {0}UPi„\ {ro} and a G {0}UPout\ {so}- The domain 
is then the set of all pairs (V, X) such that 1 — \l srX is invertible in B {^e)- In, 
the special case of unit gain we write {T) — (T, 1). For fixed V, the map 
Te (V,-) is a non- commutative fractional linear, or Mobius, transformation. 
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Lemma 14 (Siegel) If S — ^ ^ ^ ^ be a unitary operator on the di- 
rect sum Sji © Sj2 of two Hilbert spaces with \\A\\ < 1. Let $t {X) = D + 
CX (1 - AXy^ B withX e Dom{<i>s) whenever \\X\\ < 1. ForX,Y E Dom{<^s) 
we have the Siegel identities. 

$s (^)^ $s - 1 = (1 -xUt)"^ (xty - 1) (1 - AF)B, 
$5 (X) $5 (r)^ - 1 = ct (1 - xA)"^ (xr^ - 1) (i - A^r) ct. 

A proof can be found in [24] , 

Corollary 15// 5=^^ i^^*'^ unitary and 1 — yl is invertible, then D + 

C {1 — A) ^ B is unitary. 

For consistency, we would hope that the reduced model matrix belongs to 
the class of model matrices with multiplicity space one dimension lower. This 
we now show to be the case. 

Lemma 16 Let V be the model matrix determined by the operators (S,L,H). 
Then the reduced model matrix V'"'^'* obtained by eliminating the edge (ro, Sq) is 
determined by the operators (S™'', \J''^'^^ y\red\^ where 

S'"'"'' — S + S fl - S S 

L'"'"^ — L + s n - S L 

H^ed = H+ ^ ImLtS,.„(l-S.„.J-iU, 

S&Vaut 

for r G Vin\ {ro} and s G Vout\ {so}- 

Proof. The identifications VJi^'^ = S™'* ~ Ssr + '^sra (1 ^ Ssq^o)^^ ^sor and 
Mlf = Lf^ = Ls + Ssro (1 - Ssoro)"^ L^o are immediate. Unitarity of S'^'^ 
follows from the above corollary to the Siegel identities. We next check that 

\/rcd — I lodtcrcd tfp^p 

Vor - l^seVa^t\{sa}^s i^^.-tiere 

El lodtciod _ I tcrod I if C;t ^"^qt Crod 

■-s -'sr / , '-s~^sr ' / , '-sq V sorg/ ~'sro'^sr 

SePout\{so} SeVaut\{so} SeVaut\{so} 

and we use the simplification 

Ect ciod _ f S''^ S + S''^ S H-S )~^S 

S£Vaut\{so} SeVaut\{so} 

'^soro^^'or + (1 ^ Sj^^^^^Ssoro ) (1 ~ SsQro) ^so'r ^ 
{,^~'toroi^ ^ ~'soro) + (1 ~ S^jj^jj Ssoro )) (1 ~ S^oro) ^sor 



(- 
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so that 

El redfcred _ | tcred , i f (-1 _ c Q 

'-S ■-'sr ~ / , '-s-^sr ^ ""so V-^ ■-'soroJ ■~>sor- 

sePout\{so} seVout\{so} 

On the other hand 

= Vo.+Vo.„(l-V,„,J-W,,„, 

S&Vaut SeVaut 

= - \^ itS™'^ - Lt S - L^^ S (1~S V^S 

— / , '-s-^sr "-So ■-'So'" '-so~'soro \^ -^soro) -^sor 
sG7'out\{so} 

— _ \ " I tcred _ I t ('1 — «; "1"^ 9 = _ \^ I redfcred 
~ '-B-^sr '-So y- ■-'soro) ■-'sor — / , -^sr ■ 

Finally we must check that V^°'' = -^J2sev,^,\{so} ^T'^^'K"'^ - Let us 

use this as the definition of H'''^'^ then we have 



seVout\{so} 

= HH-^ E Liu)- E LtS..„ (1 - S.„.J-^ U 
sePout seVoM 



\ ^ I redfi red 
"2 / y s '-S 

seT'outXiso} 



-iH 



2 '-so '-So 2 5-^ '-s^sro (1 Ssgr-o) ^sq 

SSPoutAlso} 

+ 2 23 '-So ~ ^soro) ^sro'-s — L^oSsoro (1 ~ Ssgro) 

seVont\{so} 

+ 2 $Z '-So ~ ^so^o) ^sro^sro (1 ^ Sspro) ^sq- 

seVout\{so} 

We collected together several terms to get — ^Lj^^XL^o where 



^ ~ 1 -|- 2Ss(jro (1 Ss(jro) (l ^Soro) Ssro^sro(l Sggro) 

seVontXiso} 

= 1 + 25ggrg (1 — Ssoro) ~ (■'- ~ ^soro) ~ ^soro^soro) (l ~ Sgg^o) 
— (l ~ ^soro) (Ssoro ~ ^soro) ~ ^soro) 



'soro 

It follows that 



H-<i = H + Im E "-JSsro (1 - S«oro)"' '-so + (1 - S^^.J"^ L, 
se7'out\{so} 
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Note that Im L^^ (1 - Ss„r„) ^ \-so = Im L|„Ss„ro (1 - ^soro) ^ ^sg we obtain H'^^'^ 
as the self-adjoint operator given in the statement of the lemma. ■ 

Lemma 17 Let ei = (si, ri) and 62 = (s2, be a pair of edges in a network 
then 

Te^ O = Te^ O Te^ = J^ei®e^- (14) 

Proof. Without loss of generality suppose that ei = (1, 1) and 62 = (2,2), 
then for a,/3 ^ {1,2} 

{Te, O Te,y)^p = yap + Va2 (l - V22) V2/3 

where Vq,/? = \/a0 + V„i (1 — Vn)"^ Vi^, and so we have 

(T r.T \r\ —\/„-L- (Val'Va2) 7/ V„i 



a2 



with 



I-V22 1-Vii I-V22 

1-Vn -V12 

-V21 1 - V22 

The expression is clearly symmetric under interchange of 1 and 2, and corre- 
sponds to the double edge elimination. ■ 

This implies that the order in which we apply zero time delay limits to elim- 
inate multiple internal channels docs not in fact matter, and can be combined 
simultaneously. In this manner, every quantum network may be reduced to a 
single Markovian component in a unique well-defined algebraically manner by 
eliminating the time delays in all internal channels by means of the map 

^ = °ee£int-^e : on (f), ^net) ^ ^ % ^xt) • 

If we decide to eliminate all internal channels by making the connections, 
inserting a gain matrix X, and taking the zero time-delay limit, then the re- 
sulting network will have model matrix J- (V, X) and of course has dimensions 
determined by the remaining (external) channels. The functorial property is 
that we are able to eliminate blocks of channels in one step, giving the same 
answer as if we performed the eliminations one-by-one. 

5 Physical Applications 

We shall be interested in the situation where we eliminate all the internal chan- 
nels (total multiplicity space .^i) leaving only the external channels (total mul- 
tiplicity space With respect to the decomposition R = .ftint®-^xt, we write 
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the operators (S, L, H) of the network as 



/ Si, 






> u ) 


I Sei 









The feedback reduced model matrix may be conveniently expressed as 
^ (V, = V„;3 + V„i {rj - Vii)-^ Vi^ 

for a,/3 € {0, e}, where rj is the (unitary) adjacency m,a,trix 



Vsr 



1, if (s,r) is an internal channel, 
0, otherwise. 



Here the "gain" tj is the set of instructions as to which internal output port 
gets connected up to which internal input port. We, of course, have = 1 if we 
match up the labels of the input and output ports according to the connections, 
however, it is computationally easier to work with a general labelling and just 
specify the adjacency matrix. The reduced model matrix V''^ obtained by elim- 
inating all the internal channels is determined by the operators (S'''*^, L'''"^, H'''''^) 
given by 



gred 
l^red 
|_jred 



— See "1^ Sei Sii) Sie, 

= Le + Sei (?7 - Sii) ^ Li, 

= H+ ^ ImL]S,i(r?-Sii)-'Li. 



5.1 Quantum Systems in Feedforward 

Our first application is to derive the formula for systems in series. This is the 
simplest nontrivial example of a quantum network. 



figure 3: Systems in series 

We begin by taking (Si, Li, Hi) and (S2, L2, H2) to be the operators of the 
first and second system when considered as separate systems driven by inde- 
pendent noises. The model matrix for the network is then 

-E,-=i,2(5LjLi+^H,-) -LlSi -LtS2 
V = I Li Si 

L2 S2 
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with respect to the labelHng s = {0, 81,82} for the outputs (rows) and r = 
{0, ri,r2} for the inputs (columns). We wish to reform a feedback reduction 
wherein we connect the systems via edge e = (si,r2) and take the zero time- 
delay limit along the edge. The resulting model should then be Markovian and 
its model matrix is given by 

V —TV 

* series — e * 



E.^iAh^h+M,) -L\Si \ ^ r -452 A (Li,Si) 
L2 / \ S2 / 

E,=i,2(5L]L, + zH,) - LtS2Li -LiSi - 4S2Si \ 

L2 + S2L1 S2S1 y' 

that is, the operators determining the reduced system are 

S2S1, 
L2 + S2L1, 

Hi + H2+Im{4S2Li}. 

The Evans-Hudson maps associated with the feedforward system can be related 
to those of the individual systems via the identity 

(•) = Ci'l (•) + (M«)t42) (.) mW. (15) 

The case most familiar to the quantum optics community is two cavity sys- 
tems in cascade. Here Si =82 = 1 and Li = y^ai (i = 1,2) so we obtain 

i= i/77ai -t- ^/7^a2, H = Hi + H2 + ^^/lil2 (o-l^'i ~ aia2^ ■ This agrees with 
the calculations of Gardiner [11] for cascaded oscillators. Gardiner's derivation 
pn of the cascade rule is different from ours, but would extend to the cover the 
gauge case due to ([T5]). 




^scries 
Lfinrips 



5.1.1 The Series Product 

The rule for determining the form of the model for systems in series has been 
previously given for the Ito generator matrices where it was called the series 
product. It is related to the general question of how to "add" stochastic deriva- 
tions in order to obtain a stochastic derivation [T]. We recall its definition and 
establish its basic properties. 

Definition 18 Let Gi and G2 be ttie ltd generator matrices with the same 
multiplicity spaces, then the series product is defined to be 

G2 <]Gi ^Gi+G2 + G2nGi. 
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Lemma 19 The series product of two ltd matrices G = G2 < Gi is again 
an ltd matrix and if the Gi have parameters (S,;, L^, H) then G has parameters 
(SseriesT^-series, ^series)- If is the Galilean matrix associated with Gi then 
the Galilean matrix associated with G = G2 < Gi is M = M2M1. The series 
product is not symmetric, but is associative. 

Proof. One readily checks that the series product of two Ito matrices sat- 
isfies the conditions to be an Ito generating matrix. The specific form of the 
parameters is found by inspection. The associated Gahlean transformation is 

M = l + n(Gi + G2 + G2nGi) 

= 1 + nGi + nG2 + nG2nGi 
= (i + nG2)(i + nGi) 
= M2M1. 

To prove associativity, let us construct the (1 + n + l)-square matrix 

-LtS -iLtL 

V = I S L I (16) 

from the model parameters (S, L, H). Then the series product corresponds to 
the ordinary matrix product Vsoiics — V2V2 which is clearly associative. ■ 

Associativity means that we can extend the result immediately to several 
systems cascaded in series. The easiest way to calculate the model matrix for 
several components in series is then by the ordinary matrix product of the 
augmented matrices 

V,erie«=V„---V2Vi. 

In the lemma, the matrix V is of the type introduced by Belavkin IFI to efficiently 
capture the Ito correction as an ordinary product. Let ^ = 1 then 

VI 0/ 

we obtain an involution X* = (H^C on the space of (1 + n + l)-dimensional 
matrices. Then we are considering precisely the class of ik-unitary matrices of 
/I B A\ 

the form V = D C , that is V*V = W* = 1. In particular, the 

V 0. ^. V . 

product of two ★-unitaries is again a Tkr-unitary. 

Remark 20 Finally let us make the important remark that nowhere did we 
assume that the entries of Gi and G2 had to commute. This means that the 
series product can describe not only forward from one component system to an 
independent system, but also feedback into itself as well. This is captured in the 
picture below. 
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5.2 Beam Splitters 

A simple beam splitter is a device performing physical superposition of two 
input fields. It is described by a fixed unitary operator T=f^ ^ j € U {2): 



A2 



a (3 



Ai 
A2 



This is a canonical transformation and the output fields satisfy the same canon- 
ical commutation relations as the inputs. The action of the beam splitter is 
depicted in the figure below. On the left we have a traditional view of the two 
inputs (^1,^2) being split into two output fields {Ai,A2). On the right we have 
our view of the beam splitter as being a component with two input ports and 
two output ports: we have sketched some internal detail to emphasize how the 
scattering (superimposing) of inputs however we shall usually just draw this as 
a "black box" component in the following. 





Figure 4: Beam-splitter component. 

Our aim is to describe the effective Markov model for the feedback device 
sketched below where the feedback is achieved by means of a beam splitter. 
Here we have a component system, called the plant, in-loop and we assume that 
it is described by the parameters (So, Lq, Hq). Markovianity here corresponds 
to the limit of instantaneous feedback. 



lout 





plant 








out 


■0 





Figure 5: Feedback using a beam-splitter. 
It is more convenient to view this as the network sketched below. 
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-0^*3 S3o- 



plant 



-0S2 r-20- 



beam splitter 



Figure 6: Network representation. 

Here we have the pair of internal edges (s2,r3) and (53, r2). The model 
matrix for the network is 



V = 



V 



mio-iHo 













Til 


T12 








T21 


^22 





Lo 








So 



with respect to the labels (0,si, 82,53) for the rows and (0,ri,r2,r3) for the 
columns. Here we have 



Sii 



Li 



T22 
So 

Sei = (Ti2,0), 

Lo 




S. = 



T21 



See = Til, 



0, r? 



1 

1 



Note that the adjacency matrix has row indices (52,53) and columns indices 
(?'2,?'3) labelling the internal ports, and that the edges are the off diagonals 
(52, rs) and (53, r3). We may perform the two eliminations simultaneously to 
obtain 



Sped 



Lred — 



ired 



T11 + 



(ri2 0) / -T22 1 

\ I -So 

Til + T12 (Sq ^ — 722) 721, 

(T12 0) / -T22 1 
V 1 -So 

T12 (1 — So722) ^ Lo, 

Ho + lm (0 Lt) /-T22 1 



T21 





Lo 



, 1 -So 
Ho+ImLj(l-Sor22)"'Lo. 




Lo 
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5.3 The Redheffer Star Product 

An important feedback arrangement is shown in the figure below. 



V-OUt 

"3 



.1 



-or2 S20- 



-OS3 r-3o 



-OS4 r4CH. 



uout 

"2 



B 



Figure 7 Composite System 



We shall now derive this system taking component A to be described 



cA 

cA 
~>2l 



'-2 



, Ha and B by 



-"33 
->43 



■-3 
■-4 



Hfi. The operators of systems 



A are assumed to commute with those of B. We have two internal channels to 
eliminate which we can do in sequence, or simultaneously. We shall do the 
latter, here we have 



Se 
Si 



S44 



21 







B 
34 



■} ^ei 

S- = 



cA 
~>12 



cA 
->22 







->43 





and 



Li 



The reduced operators are therefore 



■-2 
■-3 



1 

1 







sf^ 



cA 



+ 



cA 
■^12 



-543 



cA cB 
-^12-^33 



cA cB \ 
-^22-^33 J 



cA 
•^22 


1 ^ 


)-'( 


' cA 
-^21 




1 


~^33 , 







Sf4 ) 



cB /'1 _ cA cB \-l cA 
-'43 \^ -522 •^SSj 



cA 
•->2\ 



cA 
-J 12 



cA cBV 
■^22-^33) 



-5 34 
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S44 + (1 — 5^2^^) 5^2^^ 
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y Lf + S4I (1 — S^jSfa) L^ + SgSt, (l — S^jSfa) Lf y 

H. = + Hs + Im {if (1 - SiSf2)"' + (l - SfaS^)"' S^L^ 

+'-2 V-*^ ^ S^2S£3) Lf + '-2 ^ (1 — S^2Sfi) 

+'-1 5^2 (1 — SgStz) Lf + S^2 (1 ^ SfjS^2) SfjL^ 

I I Sfc-S _ CA cS ^-l cA I B , I BtcS ('1 _ cA cS ^-l 1 a\ 
+ ■-4 -^43 \^ -'22-'33j -^22 "-3 + "-4 -^43 l-"^ -'22-'33j "-2 t ■ 



6 Topological Rules for Quantum Networks 

We wish to state algebraie rules for construeting the operators Snot, Lnct, H„ct 
covering all the examples of feedbaek networks studied so far. For a given 
network, we denote by Vout the set of output ports and V^^l the subset of 
output ports having external output. With each i € "Pout there is an associated 
port operator Li. Similarly we have Pin and P^f* for the inputs. 
For 7 = an ordered pair of ports, we set 

if j e Vout and i G Vin and there is a feedback connection 
from j to i; 

if i G Vout and j G Vin and both i and j are in the same plant 
and where S,j is the scattering component from j to i; 

otherwise. 

More generally, if 7 is an ordered sequence (paths) of ports, alternating between 
input and output ports, we define inductively. If 7 is the concatenation of 
7i followed by 72 then 

Finally, we set F to be the set of all paths going from j to i. Then for 

p -pcxt • ri -next 

(Snet)jj = X]7er(i,i) ^7' 

(Lnet)j = J2keVoM ^'rer(i,k) Sj^k', (17) 
Hnet = Ho + J2i,keVout ^7er{i,k) Lj^S^Lfe, 



S-v = < S 



1, 



where Hq is the sum of the individual component systems Hamiltonians. 
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